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Sh 1 ABSTRACT 



We have used the 2-degree Field (2dF) instrument on the Anglo-Australian Telescope to ob- 
tain redshifts of a sample of z < 3, 18.0 < g < 21 .85 quasars selected from Sloan Digital Sky 
Survey (SDSS) imaging. These data are part of a larger joint programme between the SDSS 
and 2dF communities to obtain spectra of faint quasars and luminous red galaxies, namely 
the 2dF-SDSS LRG and QSO Survey (2SLAQ). We describe the quasar selection algorithm 
and present the resulting number counts and luminosity function of 5645 quasars in 105.7 
deg 2 . The bright end number counts and luminosity function agree well with determinations 
from the 2dF QSO Redshift Survey (2QZ) data to g ~ 20.2. However, at the faint end the 
2SLAQ number counts and luminosity function are steeper (i.e. require more faint quasars) 
than the final 2QZ results from Croom et al. (2004), but are consistent with the preliminary 
2QZ results from Boyle et al. (2000). Using the functional form adopted for the 2QZ analysis 
(a double-power law with pure luminosity evolution characterized by a 2nd order polynomial 
in redshift), we find a faint end slope of f3 = —1.78 ± 0.03 if we allow all of the parameters 
to vary and (3 = —1.45 ± 0.03 if we allow only the faint end slope and normalization to 
vary (holding all other parameters equal to the final 2QZ values). Over the magnitude range 
covered by the 2SLAQ survey, our maximum likelihood fit to the data yields 32 per cent more 
quasars than the final 2QZ parameterization, but is not inconsistent with other g > 21 deep 
surveys for quasars. The 2SLAQ data exhibit no well defined "break" in the number counts or 
luminosity function, but do clearly flatten with increasing magnitude. Finally, we find that the 
shape of the quasar luminosity function derived from 2SLAQ is in good agreement with that 
derived from type I quasars found in hard X-ray surveys. 

© 2005 RAS Key words: quasars: general - galaxies:active - galaxies:Seyfert - cosmology:observations 

- surveys 



2 G. T. Richards et al. 



1 INTRODUCTION 

We have merged the high-quality digital imaging of the Sloan Dig- 
ital Sky Survey (SDSS; York et al. 2000) and the powerful spectro- 
scopic capabilities of the 2-degree Field (2dF) instrument (Lewis 
et al. 2002) to conduct a deep wide-field spectroscopic survey for 
quasars and luminous red galaxies (LRGs), i.e. the 2dF-SDSS LRG 
and QSO Survey (hereafter referred to as "2SLAQ"). The combi- 
nation of these facilities allows us to probe substantially deeper 
than either the SDSS or 2dF surveys can individually. This paper 
describes the first results of the quasar aspect of the survey; see 
Cannon et al. (2003), Padmanabhan et al. (2004) and Cannon et al. 
(2005) for a discussion of the LRG component of the survey. 

The 2dF QSO Redshift survey (2QZ; Boyle et al. 2000; Croom 
et al. 2004) was restricted to b j < 20.85. We use the SDSS imaging 
data as the input for a new survey, allowing us to probe to g — 
21.85 with typical photometric errors at the flux limit of only 7 
per cent - considerably fainter than either the i = 19.1 flux limit 
of the SDSS Quasar Survey (Richards et al. 2002; Schneider et al. 
2003) or the bj = 20.85 flux limits of 2QZ. By allocating 200 
fibres per 2dF plate to this new quasar survey (with an additional 
200 fibres going to LRGs), and extending the exposure time to 4 
hours (compared to ~ 1 hour for SDSS and 2QZ), we hope to 
obtain spectra of 10,000 quasars to g — 21.85 in the next few 
years. This paper reports on the first three semesters of data (with 
5645 quasars) and presents the z < 2.1 quasar luminosity function 
(QLF) to fainter luminosities at each redshift than either the SDSS 
or 2QZ surveys alone. 

The first determination of the luminosity function of quasars 
was by Schmidt (1968). Subsequent pioneering work was carried 
out by many groups including Schmidt & Green (1983), Koo & 
Kron (1988), Hewett, Foltz & Chaffee (1993) and especially Boyle 
and collaborators (e.g. Boyle, Shanks & Peterson 1988), with ex- 
tensions to high-z (z > 3) being provided by Warren, Hewett 
& Osmer (1994), Schmidt, Schneider & Gunn (1995), Kennefick, 
Djorgovski & de Carvalho (1995) and Fan et al. (2001b). The 
largest samples analysed to date come from the 2dF QSO Redshift 
Survey (Boyle et al. 2000; Croom et al. 2004) with 23,338 quasars. 

With the exception of variability-selected samples (e.g. 
Hawkins & Veron 1995), early QLF determinations were generally 
characterized by a strong, distinct "break" whose redshift evolu- 
tion has been the subject of much discussion. However, more re- 
cent determinations (e.g. COMBO-17; Wolf et al. 2003), while still 
exhibiting distinct curvature in a log-log plot, show less of a break 
at a specific luminosity. 

Recently, optical surveys have been supplemented by X-ray 
surveys (both soft and hard) that, when correcting for selection dif- 
ferences, can largely reproduce the optical type I QLF (e.g. Ueda 
et al. 2003; Barger et al. 2005). These X-ray QLFs have also been 
shown to exhibit a break, but generally at luminosities much fainter 
than found by optical surveys; this result suggests incompleteness 
at the faint end of optical surveys. 

As we shall see, our data are in good agreement with recent 
results for faint quasars from both the optical (e.g. Wolf et al. 2003) 
and X-ray (e.g. Barger et al. 2005). We probe nearly 1 magnitude 
deeper than 2QZ, and find that the faint end slope of the QLF is 
steeper than that of the most recent 2QZ determination and lacks a 
strong characteristic break feature, but is still better characterized 
by a double-power law than a single power-law. 

Section 2 presents a description of the imaging data and the 
sample selection. In Section 3 we describe the observations and 
data reduction. Section 4 presents the completeness corrections 



leading to the QLF presented in in § 5. Finally, § 6 presents a 
discussion of the ramifications of our work and summarizes our 
results. Throughout this paper we use a cosmology with H — 
/i 7 o70kms _1 Mpc~\ Q m = 0.3, Qa = 0.7 (e.g. Spergel et al. 
2003). 



2 THE IMAGING DATA AND SAMPLE SELECTION 

2. 1 The SDSS imaging data 

The photometric measurements used as the basis for our catalogue 
are drawn from SDSS imaging data (DR1 reductions; Stoughton 
et al. 2002; Abazajian et al. 2003), which will eventually cover 
roughly 10,000 deg 2 of sky in five photometric passbands (ugriz) 
using a large-format charge-coupled device (CCD) camera (Gunn 
et al. 1998). The photometric system and its characterization are 
discussed by Fukugita et al. (1996), Hogg et al. (2001), Smith 
et al. (2002) and Stoughton et al. (2002); the spectroscopic tiling 
algorithm is described by Blanton et al. (2003). Except where oth- 
erwise stated, all SDSS magnitudes discussed herein are "asinh" 
point-spread-function (PSF) magnitudes (Lupton, Gunn & Szalay 
1999) on an AB magnitude system (Oke & Gunn 1983) that have 
been dereddened for Galactic extinction according to the model of 
Schlegel, Finkbeiner & Davis (1998). The astrometric accuracy of 
the SDSS imaging data is better than 100 mas per coordinate rms 
(Pier et al. 2003). The SDSS Quasar Survey (Richards et al. 2002; 
Schneider et al. 2003, 2005) extends to i = 19.1 for z < 3 and 
i = 20.2 for z > 3, whereas our work herein explores the z < 3 
regime to g = 21.85 (i ~ 21.63). 

2.2 Preliminary sample restrictions 

Our quasar candidate sample was drawn from 10 SDSS imaging 
runs (see § 2.4) after having first been vetted of objects that have 
cosmetic defects (e.g. bad columns) that might cause the photome- 
try to be inaccurate. Specifically, we rejected any objects that met 
the "fatal" or "non-fatal" error definitions as described by Richards 
et al. (2002). 

We next imposed limits on the i-band PSF magnitude and its 
estimated la error of i < 22.0 and Oi < 0.2. Further magnitude 
cuts are done in the p-band (to facilitate comparisons with previ- 
ous 2QZ results in fej); the i-band cuts are primarily to reduce the 
number of objects that we have to examine initially. We also placed 
restrictions on the errors in each of the other four bands, specifi- 
cally, a u < 0.4, o g < 0.13, o r < 0.13 and a z < 0.6. These 
restrictions are designed to ensure that the errors on the magni- 
tudes are reasonably small (and thus that the resulting colours are 
accurate), but also are sufficiently relaxed that, when coupled with 
the magnitude cut in i and g, objects with quasar-like continua are 
not rejected. This tolerance is necessary since, as we go fainter, re- 
strictions on magnitude errors are effectively cuts in magnitude and 
any two such restrictions are effectively colour cuts. Note that this 
selection of error constraints effectively limits the redshift to less 
than 3, as the Lya forest suppresses the u flux at higher redshifts. 

2.3 Colour cuts 

Based on spectroscopic identifications from SDSS and 2QZ of this 
initial set of objects, we implement additional colour cuts that are 
designed to efficiently select faint quasars while maintaining a high 
degree of completeness to known UV-excess broad-line quasars. 
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An analysis of the completeness of the selection algorithm is given 
as a function of redshift and magnitude in § 14.21 

We first impose colour restrictions that are designed to re- 
ject hot white dwarfs. These cuts are made regardless of magni- 
tude. Specifically, we rejected objects that satisfy the condition: 
A&&i((B&zLC&z&zD)\ \E), where the letters refer to the cuts: 



(1) 



This is similar to the white dwarf cut applied by Richards et al. 
(2002, Equation 2) except for the added cut with respect to the g — i 
colour. 

As the targets become fainter and the magnitude errors in- 
crease, we find that maximizing our completeness and efficiency 
is best served by separate handling of bright and faint objects. The 
bright sample is restricted to 18.0 < g < 21.15 and is designed to 
allow for overlap with previous SDSS and 2dF spectroscopic ob- 
servations. The faint sample has 21.15 ^ g < 21.85 and probes 
roughly one magnitude deeper than 2QZ. These cuts are made in g 
rather than i (as the SDSS quasar survey does) since we are con- 
centrating on UV-excess quasars and would like to facilitate com- 
parison with the results from the bj-based 2QZ. The combination 
of the g < 21.85 and i < 22.0 cut will exclude objects bluer than 
a v — +0.3 (/„ oc v a ) \ however, objects this blue are exceedingly 
rare (> 3a deviations). 

Further cuts are made as a function of colour and morphol- 
ogy in each of the bright/faint samples. In general, we would pre- 
fer not to make a cut on morphology since we do not want to ex- 
clude low-z quasars and because our selection extends beyond the 
magnitude limits at which the SDSS's star/galaxy separation breaks 
down. However, Scranton et al. (2002) have developed a Bayesian 
star-galaxy classifier that is robust to r ~ 22. As a result, in addi- 
tion to straight colour-cuts, we also apply some colour restrictions 
on objects with high r-band galaxy probability (referred to below as 
"galprob") according to Scranton et al. (2002) in an attempt to re- 
move contamination from narrow emission line galaxies (NELGs; 
i.e. blue star-forming galaxies) from our target list. 

Bright sample objects are those with 18.0 < g < 21.15 and 
that meet the following conditions 
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Figure 1. Colours of spectroscopically identified 2SLAQ targets. Con- 
firmed quasars are shown in black, non-quasars in red. For reference the 
colour distribution of z < 3 SDSS-DR1 quasars (Schneider et al. 2003) are 
shown as faint grey contours/points. The dashed and solid dark blue lines 
show the u — g and g — i colour cuts for the bright and faint samples, re- 
spectively. The dotted cyan line shows the boundary of the g — i cut used 
to reject white dwarfs. 
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(3) 



in the combination A&cMB&iMC&cMD&cME, where cut A se- 
lects UVX objects, cuts B, C and D eliminate faint F-stars whose 
metallicity and errors push them blueward into the quasar regime, 
and cut E removes NELGs. These faint cuts are more restrictive 
than the bright cuts to avoid significant contamination from main 
sequence stars that will enter the sample as a result of larger errors 
at fainter magnitudes. 

FigureQshows the u — g vs. g — i colour distribution of objects 
satisfying these criteria for which we obtained new spectra. Objects 
confirmed to be quasars are shown in black, while those that are not 
quasars (mostly stars and NELGs) are shown in red. The locus of 
z < 3 quasars from SDSS-DR1 (Schneider et al. 2003) is given by 
grey contours and points. Solid blue, dashed blue and dotted cyan 
lines show the faint sample colour cuts, the bright sample colour 
cuts and the white dwarf cut, respectively. 



in the combination A&iMB&z&zlC&iMD&zME, where cut A se- 
lects UVX objects, cuts B and C eliminate faint F-stars whose 
metallicity and errors push them blueward into the quasar regime, 
and cuts D and E remove NELGs that appear extended in the r 
band. Among the bright sample objects, those with g > 20.5 were 
given priority in terms of fibre assignment. 

Faint sample objects are those with 21.15 ^ g < 21.85 and 



2.4 Sky location of imaging data 

2.4.1 2003 A and 2004A 

For the first semester both of 2003 and 2004, we used the SDSS 
imaging data (rerun 20; Stoughton et al. 2002; Abazajian et al. 
2003) in the SDSS northern equatorial scan (stripe 10) from runs 
752, 756, 1239 and 2141; see York et al. (2000) and Stoughton et al. 
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Figure 2. RA and Dec distribution of 2dF targets selected from SDSS imag- 
ing data (black) and previous 2dF observations (red). 2SLAQ observations 
are given by cyan points within 41 plates (green circles). The top two panels 
are semester A targets; the bottom two panels are semester B targets. Note 
the distortion of the coordinate system; the spectroscopic plates are actually 
circular. 

(2002) for the definition of the relevant SDSS technical terms. Run 
756 was used for the northern part of the stripe, while a combi- 
nation of the other three runs was used for the southern part of 
the stripe in an attempt to use the best quality data (typically that 
with the best seeing). The area of sky sampled was further limited 
to regions where the SDSS image quality was deemed to be good 
enough to use for targeting faint objects for spectroscopy, specifi- 
cally seeing $C l'/8 and r-band Galactic extinction ^0.2 (Schlegel 
et al. 1998). We also excluded any objects from SDSS camera col- 
umn 6, since 2dF cannot cover the full 2.5 degree wide SDSS stripe 
and column 6 has the lowest quality data of all the columns as a 
result of (relatively) poorer image quality at the edge of the cam- 
era in these early SDSS data. The final RA ranges were 123° < 
aj2ooo < 144°, 150° < ajaooo < 168°, 185° < Qjaooo < 193°, 
197° < aj2ooo < 214° and 218° < a J2 ooo < 230°. Whenever 
possible, we tried to overlap areas with existing 2QZ spectroscopy 
to limit the number of objects with b.j < 20.85 that needed spec- 
troscopic confirmation. SDSS spectroscopy limits the need for new 
i < 19.1 spectra. The two top panels of Figure|2|illustrate the area 
covered by our Semester A targets (—1.259° ^ <5j20oo ^ 0.840°). 

2.4.2 2003B 

For the second semester, our samples were limited to the following 
combination of data (run, rerun, strip, aj2uoo range): (2659, 40, 
82N, 309.20° < a.72000 < 320.34°), (2662, 40, 82N, 320.34° < 
aj2ooo < 15.08°), (2738, 40, 82N, 15.08° < a. J2 ooo < 59.70°), 
(2583, 40, 82S, 309.20° < a J20 oo < 341.08°), (3388, 40, 82S, 
341.08° < QJ2000 < 345.44°), (3325, 41, 82S, 345.44° < 
QJ2000 < 59.70°). These reruns (40 and 41) represent post-DRl 
data processing, which includes a newer version of the photomet- 
ric pipeline and improved photometric calibration. Again, camera 
column 6 was excluded and these are all equatorial scans. Note that 



there are no 2QZ observations in this range. The two bottom panels 
of Figure|2|illustrate the area covered by our Semester B targets. 

2.4.3 Sky area 

The area of sky covered by our catalogue of targets for 2003/4A 
was 159.4 deg 2 with 20228 targets and for 2003B it was 230.2 deg 2 
with 33160 total targets. Thus we have a total area of 389.6 deg 2 
and 53388 targets. Of this area, this paper concentrates on only 
those regions where we have obtained new spectra (see Figure [2}. 
In semester 2003/4A, 34 plates were observed, covering an area of 
80.82 deg 2 - as determined by the fraction of targets within the 
plate areas (11075 of 53388). In the second semester, seven plates 
were observed, covering an area of 24.9 deg 2 (3407 of 53388 tar- 
gets within the new plate areas). Note that the plates overlapped in 
2003/4A, but not in 2003B. The theoretical area for 2003B given a 
plate radius of 1.05 deg is 24.3 square degrees, which compared to 
the area estimated by fraction of targets (24.9 deg 2 ) suggests that 
our estimate of the area has a roughly 2.5 per cent error. Thus the 
area covered by new plate observations is 105.7±2.6 deg 2 . Within 
these plate centers there are 14482 targets, of which 9120 have 
spectroscopic identifications, and among those are 5645 quasars. 



3 SPECTROSCOPIC OBSERVATIONS 

3.1 The 2dF facility 

The Two Degree Field (2dF; Lewis et al. 2002) facility at the 
Anglo- Australian Telescope is a fibre fed multi-object spectrograph 
and robotic fibre positioner. The fibres are 140/im in diameter, 
which is roughly 2'.' 16 at the centre of the plate and l'.'99 at the 
edges. Two independent spectrographs use Tektronix 1024x1024 
CCDs with a range of diffraction gratings offering resolutions be- 
tween 10A and 2.2A over the optical wavelength range. During 
standard operation, 400 fibres are available for simultaneous ob- 
servation (200 per spectrograph) over a 2 degree diameter field of 
view. The system is equipped with an atmospheric dispersion com- 
pensator which enables 2dF observations to be taken over a wide 
wavelength range, by ensuring that all wavelengths from the UV to 
NIR enter the fibres. However, differential spatial atmospheric re- 
fraction distorts the field geometry and limits observations of equa- 
torial fields to ±1 hour on either side of transit. 

3.2 2dF field configurations 

The 2SLAQ survey regions are centred close to the equator and 
are 2 degrees wide in declination. To achieve optimal sky cover- 
age while still retaining a largely contiguous area, the 2dF field 
centres are placed along the central declination of the two strips; 
<5j2ooo = -00° 12'35" for the North Galactic Pole (NGP) strip and 
£12000 = -00°15'00" for the South Galactic Pole strip; see 
Each field centre is separated by 1.2 degrees, although some early 
observations of the NGP at the start of 2003 had field spacings of 1 
degree. 

The target list generated from the process described in §|2]is 
then merged with a target list of LRGs selected from the same pho- 
tometric data set (Cannon et al. 2005). The sub-samples within this 
combined data set are assigned different priorities which determine 
the likelihood of a fibre being assigned to them in the 2dF config- 
uration process. The priority values given to each sample are listed 
in Table Q where 9 is the highest and 1 the lowest priority. All 
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Table 1. 2dF configuration priorities. 



Sample 


Priority 


Guide stars 


9 


LRG (main) random 


8 


LRG (main) remainder 


7 


QSOs (g > 20.5) random 


6 


QSOs (g > 20.5) remainder 


5 


LRG(extras)+hi-z QSOs 


4 


QSOs (g < 20.5) 


3 


previously observed 


1 



available high priority targets are allocated before moving to the 
next priority level. For source densities much greater than 400 per 
2dF field, the 2dF configuration algorithm will tend to give a non- 
uniform distribution of fibres allocated to objects (Cannon et al. 
2005). Therefore the main samples of each of the LRG and QSO 
data sets were randomly sampled to a surface density of 200 per 
2dF field and these given priority 8 and 6 for the LRGs and QSOs 
respectively. The remaining sources from these main samples were 
given lower priority (7 and 5 respectively). Other sub-samples, such 
as bright QSOs and high-redshift candidates, were given lower pri- 
ority. 

For the QSO sample we used the low resolution 300B grating 
(as used for the 2QZ), but the LRG observations required the use of 
the higher resolution 600V grating. Therefore, one of the 2dF spec- 
trographs is configured with a 300B grating (spectrograph 1) while 
the second (spectrograph 2) is configured with the 600V grating. 
On each 2dF field plate of 400 fibres, each block of 10 fibres (a 
retractor block) goes to an alternate spectrograph, so that 200 fibres 
on each plate are available for the QSOs and 200 for the LRGs. 
2dF fibres are limited to a maximum off-radial angle of 14 degrees, 
therefore there are 20 small triangles surrounding the edge of the 
2dF fields that are inaccessible to the QSO spectrograph covering a 
total area of 0.43 deg 2 . The angular completeness function defined 
by this complex field pattern is not relevant to the QLF analysis be- 
low, but it is critical to accurate measurements of clustering. Of the 
200 fibres available for the QSOs, 20 were allocated to positions 
known to be blank sky to be used for sky subtraction. 



3.3 2dF observations and data processing 

Observations started in March 2003. Each 2dF field was observed 
for a minimum of four hours (more if weather was poor). These 
four hours were split over two nights to minimize the effects of 
changing atmospheric refraction. The 300B grating used gives a 
dispersion of 4. 3 A pixel -1 and an instrumental resolution of 9A. 
The spectra cover the range 3700-7900A. The data were reduced 
in real time using the standard 2dFdr pipeline (Bailey et al. 2003, 
MNRAS submitted). The exposure times increased if the condi- 
tions meant that a pre-defined completeness limit (80 per cent) was 
not met. Any source that has a high S/N spectrum and a high- 
confidence identification after the first night of observation has its 
fibre assigned to previously unallocated sources for future observa- 
tions of the field. 

The identification of QSOs and measurement of redshifts was 
done using the AUTOZ code that was developed for the 2QZ (see 
Croom et al. (2001, § 3.1) and Croom et al. (2004, § 2.3.1) for 
details). All spectra are then checked by eye to confirm the iden- 
tifications. Since spectroscopic processing is the same as that used 
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Figure 3. Spectroscopic (dashed line) and coverage (solid line) complete- 
ness fractions as a function of magnitude in 0.25 mag bins. These correc- 
tions were applied (in conjunction with the photometric completeness cor- 
rection [Fig.|4]) to determine the corrected number counts. The dotted verti- 
cal lines show the boundaries of the SDSS (i = 19.1, g ~ 19.32) and 2QZ 
surveys (bj ~ g = 20.85). The discontinuity in coverage completeness at 
g ~ 20.5 is caused by a prioritization of targets fainter than this limit; see 
Table 1. 



for 2QZ spectra (e.g., quasars must have broad [> 1000 km s -1 ] 
emission lines), we treat 2SLAQ selected objects with 2QZ spectra 
as if they were observed as part of the 2SLAQ programme. 



4 COMPLETENESS CORRECTIONS 

In this section we explore and quantify the various effects that will 
bias the quasar number counts and luminosity function. In partic- 
ular, we address the photometric, coverage, spectroscopic and cos- 
metic defect incompleteness of our sample. In addition, we investi- 
gate the difference between g and bj magnitudes, Eddington bias, 
morphology bias and the effects of variability. 



4.1 Coverage and spectroscopic completeness 

We have not obtained spectra of all our quasar candidates in the 
105.7 deg 2 analyzed in this paper. Thus we must compute the "cov- 
erage" completeness of our sample, which multiplied by the area 
yields the effective area of the survey. Since we are combining data 
from three distinct surveys (SDSS, 2QZ, and 2SLAQ) in order to 
increase our dynamic range, it is necessary to compute this cor- 
rection as a function of magnitude. The coverage completeness is 
computed under the assumption that the fraction of objects that re- 
main unobserved (at a given magnitude) will be quasars at the same 
rate as those that are observed. This assumption is reasonable given 
that the objects observed are chosen at random. Figure|3|shows the 
coverage completeness (solid line) that we compute as a function 
of magnitude. 

In addition to the coverage completeness, we must correct 
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for those cases in which our spectroscopy does not yield an un- 
ambiguous identification. Assuming that the fraction of unidenti- 
fied objects will be quasars at the same rate as those among iden- 
tified objects (as a function of magnitude), we derive a spectro- 
scopic incompleteness as shown by the dashed line in Figure [3] 
This assumption arguably may tend to overestimate the number of 
unidentified quasars since the spectroscopic completeness may ad- 
ditionally be a function of redshift (because of emission line ef- 
fects which generally facilitate quasar identification) and that any 
completeness determination is surely to be a lower-limit. However, 
our spectroscopic completeness is generally high (70 per cent at 
the faint limit, 90 per cent brighter), thus any second-order correc- 
tions will have a minimal impact. Furthermore, comparison with 
supplementary identifications based solely on photometry and pho- 
tometric redshifts (§[5) suggests that this assumption is reasonable. 
In practice we have treated the spectroscopic completeness as if the 
unidentified objects simply had not been observed, which facilitates 
the application of these corrections to our model of the luminosity 
function. 
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4.2 Photometric completeness 

The incompleteness of our sample due to colour cuts is a strong 
function of both redshift and magnitude since the colours of quasars 
change significantly with redshift and fainter quasars have larger 
errors. We quantify this incompleteness by running our selection 
algorithm on a sample of simulated quasars that were designed to 
test the SDSS's quasar target selection algorithm; see Fan (1999) 
and Richards et al. (2005). The primary independent variable in the 
simulations is the spectral index distribution, which is given by a 
Gaussian distribution with a v — —0.5 ± 0.3 (/„ oc v a "), which 
is in good agreement with the composite SDSS quasar spectrum 
given by Vanden Berk et al. (2001). Blueward of the Lya emission 
line we instead use a spectral index of a v = —1.5 ± 0.17, con- 
sistent with Telfer et al. (2002); this spectral index is taken to be 
uncorrelated with the optical/UV spectral index. Only the spectral 
index, the Lya equivalent width and the Lya forest strength vary; 
all other emission lines are fixed relative to Lya. 

Figure |4] shows the selection completeness to these simulated 
quasars as a function of redshift and g magnitude. Two represen- 
tative ranges are shown, with bins 0.25 mag wide centered on 
g = 20.775 and 21.525. The g — 20.775 completeness curve 
(dashed line) is representative of the "bright" sample, whereas the 
g — 21.525 curve (solid line) is representative of the "faint" sam- 
ple (except for the faintest bin since it extends to g = 21.9 and the 
sample only goes to g — 21.85). Incomplete redshift regions occur 
when photometric errors are large and/or emission/absorption lines 
bring the colours of the quasars near/across the colour cuts (e.g. 
Richards et al. 2001). 



4.3 Correction for cosmetic defects 

Certain cosmetic defects within the imaging data cause quasars to 
be missed from our sample. Thus, we need to make a correction for 
cosmetic defects in the SDSS data, specifically for those objects 
that fail the fatal/non-fatal error tests (Richards et al. 2002). One 
way to quantify this is to assume that any cosmetic defects that pre- 
vent the selection of a particular quasar in the SDSS imaging are 
unlikely to have been present in the 2QZ imaging inputs. With the 
exception of blended objects, this assumption should be roughly 
true. Thus we match the NGP sample of quasars from the 2QZ to 



Redshift 

Figure 4. Completeness as a function of redshift and g magnitude based on 
simulated quasars. Representative magnitude ranges are shown for "bright" 
and "faint" samples with redshift intervals of 0.05. 



our initial catalogue of Semester A targets (with only the fatal and 
non-fatal errors, i < 22 and cr; < 0.2 cuts applied). Since the 2QZ 
only went to bj = 20.85, the i magnitude cut should not cause us to 
lose many quasars; however, the fatal and non-fatal error cuts (i.e. 
cosmetic defects) will cause quasars to be lost. The fraction of 2dF 
quasars that are not among our initial SDSS-imaging selected sam- 
ple gives us an estimate of the fraction of quasars that are missed 
due to cosmetic defects. We find that this fraction is ~ 5 per cent. 
A similar fraction is derived by Vanden Berk et al. (2005) based 
on an empirical analysis of the point-source completeness of the 
SDSS quasar catalogue. We apply this correction independent of 
magnitude 1 and redshift in addition to the coverage, spectroscopic 
and photometric completeness corrections described above. Losses 
due to blending of sources will increase this completeness correc- 
tion; for our purposes such losses are assumed to be smaller than 
the other corrections that we apply. 



4.4 Eddington bias 

Eddington bias is the distortion of the object number counts as a 
function of magnitude that occurs when photometric uncertainty 
causes errors in distributing sources into their proper magnitude 
bins. The relationship between the observed and actual differential 
number counts, A(m), is given by Peterson (1997): 



A(m) « A obs (m) 



log e 



(4) 



where a is the Gaussian error in the magnitude, k is defined 
by the integrated number counts relation N(m) oc C10 Km , and 
log A(m) = C + Km. If the product of the slope and the error 

1 But note that Vanden Berk et al. (2005) find that this completeness is a 
function of magnitude; however, the completeness has not been determined 
at the faint limits to which we are probing, so we assume a uniform value. 
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Figure 5. Left: Number of spectroscopically confirmed quasars classified by PHOTO as stellar (black) and as extended (grey) as a function of Bayesian galaxy 
probability (Scranton et al. 2002). Right: Fraction of point-like quasars (as determined by the Bayesian analysis) that are mis-classified by PHOTO as extended 
- as a function of magnitude. 



(an) increases with magnitude then the observed slope is steeper 
than the intrinsic slope; for decreasing an the observed slope is 
flatter than the intrinsic slope. For our sample the correction term 
in brackets above is > 0.98 for all magnitude bins, thus we have 
applied no correction for Eddington bias. 

4.5 Morphology bias 

Our sample includes objects that the SDSS photometric pipeline 
(PHOTO; Lupton et al. 2001) classifies as extended. The rationale 
for this decision is summarized in Figure[5]which shows that at the 
faintest limits of our survey, a significant fraction of point sources 
are mis-classified by the photometric pipeline as extended (assum- 
ing that the Bayesian analysis of Scranton et al. 2002 represents 
ground truth). The right-hand panel shows that this is a function of 
magnitude. The left-hand panel shows the Bayesian galaxy proba- 
bility distribution for both point-like (stellar) and extended quasars 
as classified by PHOTO. 

The inclusion of extended sources can lead to a bias. Specif- 
ically, since many of the Semester A targets have been observed 
as part of the 2QZ and since the 2QZ did not target extended 
sources, our new observations will be preferentially biased towards 
extended sources. Thus our corrections from the number of objects 
observed to the number of objects targeted may be skewed since it 
assumes that new observations will yield quasars at the same rate 
as old observations. However, we find that, although the contami- 
nation among extended sources is larger than for point sources, the 
shape of the corrections as a function of magnitude are not sig- 
nificantly different and thus our analysis of the shape of the QLF 
should not be adversely affected. 

4.6 g vs. bj 

To properly compare our 2SLAQ results to those of the 2QZ, we de- 
termine the relationship between the SDSS g band and the bj band 
used by 2dF. Figure |5| shows the two transmission curves, which 



are quite similar. The bj curve was kindly provided by Paul Hewett 
(2004, priv. comm.). The g curve is as taken from the SDSS web 
site 2 - except that it has been converted from 1.3 to 1.0 airmasses 
(to match the bj curve). In FigureQ we plot the g — bj magnitude 
difference versus bj for all of the 2QZ quasars in our sample. This 
plot shows that even considering the scatter in g — bj, the g band 
magnitude limits of our current sample completely encompass the 
2QZ quasars. 

To convert bj to g we simply compute the median g — bj dif- 
ference, which is shown as a function of bj by points connected 
by solid lines in Figure Q The median for the whole sample is 
(g — bj) = —0.045, with no significant dependence on bj. Given 
the empirical similarity of the g and bj magnitudes, and that the 
error in the computed median is of order the median itself, we have 
simply taken bj as an exact surrogate for g in our comparison of 
the number counts and luminosity functions. 

Much of the scatter between bj and g is caused by variability 
in the ^20 years between the epochs when the bj and g data were 
taken - in contrast with the simultaneity of the SDSS 5-band imag- 
ing data. The scatter in g — bj is a g -bj = 0.25 at bj — 18.475 
and Og-bj = 0.38 at bj = 20.725. At least 0.15 mag of this er- 
ror is due to photometric error in b.j (Croom et al. 2004, Fig. 9); 
roughly 0.02 and 0.035 is due to photometric error in g. Thus, 
most of the scatter (roughly a ~ 0.2) is thus caused by variabil- 
ity. Variability adds uncertainty to the magnitude distribution in the 
same manner as photometric errors and thus can modify the num- 
ber counts through Eddington bias. Proper treatment of variability 
in light of quasar number counts is complicated, ideally using long 
terms averages of the quasars under consideration. However, we 
can estimate the effect that variability has on the slope of the num- 
ber counts. If the variability amplitude is constant with magnitude, 
then variability will cause a slight flattening of the observed distri- 
bution due to the number counts being steeper at the bright end than 
at the faint end. For (j va r = 0.2, at g ~ 18.5 the number counts 

2 http://www.sdss.Org/dr3/instruments/imager/filters/g.dat 
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Figure 6. UKST bj (P. Hewett 2004, private communication) and SDSS g 
transmission curves for 1 airmass and normalized such that J S\ 4^ = 1. 
Both curves are given in terms of detector quantum efficiency, which means 
the g curve is shown as published, but the bj curve has been multiplied by 
wavelength. 
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Figure 7. g — bj magnitude differences versus bj, where g is from the 
SDSS and bj is from the 2QZ. The dotted lines show the g magnitude 
limits of the 2SLAQ sample. The dashed lines show the bj magnitude limits 
of the 2QZ survey. The points connected by a solid line give the median 
g — bj as a function of bj in 0.5 mag bins. The median over all bj is 
(g — bj) = —0.045. Neither g or bj are extinction corrected in this plot. 



will be over-estimated by ~8% and at g ~ 20.7 they will be over- 
estimated by ~1%, which produces a negligible (~2%) change in 
slope over this range. 



5 NUMBER COUNTS AND LUMINOSITY FUNCTION 

5.1 Redshift and absolute magnitude distributions 

Having discussed the various completeness corrections, we can 
now determine the number counts and luminosity function of our 
sample. Figure [8] shows the M B vs. redshift distribution of spec- 
troscopically confirmed 2SLAQ, 2QZ and SDSS quasars in our 
sample - within the boundaries of new plate observations (105.7 
deg 2 ). The absolute g magnitude, M g , is computed using lumi- 
nosity distances in the cosmology given in § 1 according to the 
prescription of Hogg (1999) and with the (albeit poor, but com- 
monly used) assumption of a universal power-law continuum of 
a v = — 0.5 (/„ oc u a ). 3 Black, cyan and blue points represent 
new 2SLAQ quasars, previously confirmed 2QZ quasars and pre- 
viously confirmed SDSS quasars, respectively. Dashed red lines at 
g = 18.0 and 21.85 demarcate the g magnitude boundaries of our 
sample. In addition we show the g ~ bj — 20.85 limit of the 2QZ 
survey. The histograms to the left and bottom of the figure show 
the one-dimensional distribution of sources in M g and redshift. We 
further overlay a grid which highlights the magnitude and redshift 
bins that were used in the construction of the Croom et al. (2004) 
QLF and will also be used for determining the binned 2SLAQ QLF. 



5.2 Number counts 

Figure [5] shows the differential number counts as a function of g 
magnitude in bins of 0.25 mag, both corrected (solid circles) and 
uncorrected (open circles) for the various sources of incomplete- 
ness (error bars are Poisson). Number counts from 2QZ (Croom 
et al. 2004) are shown in red for comparison. This diagram only 
includes quasars with M g < —22.5 and 0.4 < z < 2.1. 

From Figure|5|we see that to g ~ 20.2 the agreement between 
2SLAQ and 2QZ is quite good, but there is a discrepancy between 
the two studies at the faint end: 2SLAQ suggesting a higher density 
of faint objects than 2QZ. We note that the shape of the distribu- 
tion is clearly better fit by a double power law than a single power 
law (demonstrating the turnover in the distribution towards fainter 
quasars), but that the change in slope is more subtle than the dis- 
tinctive "break" near g ~ 19.5 that is sometimes found in such 
analyses (e.g. Boyle et al. 1987). This behaviour is qualitatively 
consistent with that found by Wolf et al. (2003) from the COMBO- 
17 survey and is inconsistent with the single power-law form found 
in variability selected samples (e.g. Hawkins & Veron 1995, but see 
Ivezic et al. 2004). 

We have shown (as open circles) the raw number counts to 
give the reader an idea of the absolute lower limits on the points and 
the size of the completeness corrections that have been applied. The 
coverage corrections are straightforward and should be fairly robust 

3 Ideally we would determine a spectral index for each individual object. 
However, this requires better spectrophotometry/photometry at the faint end 
than 2SLAQ provides. Fortunately, the errors induced by assuming a fixed 
spectral index are mitigated by the z < 2.1 nature of our analysis (the 
errors increase with redshift) and the fact that the majority of quasars have 
roughly this spectral index. 
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Redshift 

Figure 8. Absolute g magnitude versus redshift for all confirmed quasars in 
our sample. Blue crosses are SDSS quasars, cyan crosses are 2QZ quasars, 
while black dots represents quasars confirmed by 2SLAQ. The bottom and 
side histograms show the ID distributions of M g and redshift. The dashed 
red lines show the bright and faint magnitude limits of this survey (g = 
f8.0 and g = 21.85) and the approximate limit of the 2QZ survey in g 
(bj = 20.85). The grid of grey lines outline the bins used for determining 
the quasar luminosity function. 



10 



QO 



C\2 



QO 
T3 



QO 



1 1 1 1 I 1 ' ' ' 

o 2SLAQ raw 

• 2SLAQ corrected 

• 2SLAQ w/ photo-z 

• 2QZ/6QZ rjT 8 " 1 




HB-HB-j 



19 



20 



21 



22 



Figure 9. 2SLAQ number counts (black circles) compared to 2QZ number 
counts (red filled triangles) from Croom et al. (2004). The 2SLAQ number 
counts are given both as raw (i.e. observed; open circles) and corrected 
(filled circles) number counts. 2SLAQ number counts are also given (as grey 
circles) after including photometrically identified quasars (with photometric 
redshifts) from Richards et al. (2004). Quasars are restricted to those with 
M g < —22.5 and 0.4 < z < 2.1 (3889 quasars) for comparison with 
the 2QZ number counts. The dotted vertical line marks the dividing line 
between the 2SLAQ bright and faint samples. 
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Figure 10. Comparison of 2SLAQ quasar number counts to previous deep 
samples. 2SLAQ quasars are limited to M g < —23 and 0.6 < z < 2.2 
to mimic the exclusion of extended sources (which mostly have z < 0.6 or 
M g > —23). Open squares indicate points from Boyle and collaborators, 
specifically Boyle, Shanks & Peterson (1988) [cyan]; Boyle, Fong, Shanks 
& Peterson (1990) [blue]; Boyle, Jones & Shanks (1991; BSJ91) [magenta]. 
Open triangles refer to Koo & Kron (1988), where the red triangles are for 
0.9 < z < 3.0 and the grey triangles are for z < 2.2 (and z > 0.9) 
as given by Table 8 in B JS9 1 . Filled magenta pentagons refer to Marano, 
Zamorani & Zitelli (1988), as given by Table 8 in BJS91. The filled blue 
pentagon is derived from Zitelli et al. (1992). 



(perhaps less so in the 3 faintest bins due to the more restrictive se- 
lection criteria and larger photometric error). In fact, we could have 
simply corrected the effective area as a function of magnitude and 
shown the (more complete and much smoother) area-corrected raw 
counts. However, as we are splicing together three samples (SDSS, 
2QZ, and 2SLAQ) to provide spectroscopic coverage of our targets, 
it seems appropriate to fully disclose the magnitude dependence of 
the coverage completeness within the 105.7 deg 2 area covered by 
the 2SLAQ plates. As a check on our correction terms, we have 
also matched our unobserved and unidentified objects to the pho- 
tometric quasar candidate catalogue of Richards et al. (2004a), in 
attempt to "observe" a larger fraction of our quasar candidates (to 
g < 21). The objects from Richards et al. (2004a) are expected 
to be 95 per cent accurate (averaged over all magnitdues) with re- 
spect to quasar classification, with 90 per cent having photometric 
redshifts correct to \Az\ ± 0.3 for the redshift range considered 
here (Weinstein et al. 2004). The result of including these photo- 
metric identifications is shown by the grey points in Figure [5] and 
lends credence to the steeper faint-end number counts relations that 
we derive solely from our (completeness-corrected) spectroscopic 
sample. This comparison is meant purely as a sanity check. The 
differences between the spectroscopic (black cirles) and photomet- 
ric (grey circles) number counts are consistent with the expected 
decrease in efficiency of the R04 photometric catalog with fainter 
magnitude, thus supporting the accuracy of our completeness de- 
terminations (and our corrected number counts). 

We further compare our results to a number of other sam- 
ples of faint quasars that pre-date the 2QZ sample. This compar- 
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Table 2. Cumulative number count comparison. The cumulative number 
counts are shown for the Boyle et al. (2000), Croom et al. (2004) and 
2SLAQ maximum likelihood parameterizations for 16 < g < mag and 
0.3 < z < 2.2 in unit of counts per square degree. 
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ison is shown in Figure flQl Here we have restricted our sample 
to 0.6 < z < 2.2 and M g < —23 to best mimic the limits of 
these previous surveys which generally excluded extended sources 
(which typically have z < 0.6 or M g > —23). We specifically 
compare our 2SLAQ results to the samples of Boyle et al. (1988), 
Koo & Kron (1988), Marano, Zamorani & Zitelli (1988), Boyle 
et al. (1990), Boyle, Jones & Shanks (1991) and Zitelli et al. (1992), 
where Table 8 in Boyle et al. (1991) is the source of the Boyle et al. 
(1990), Zitelli et al. (1992), and Koo & Kron (1988) [z < 2.2] 
points. The redshift ranges and magnitude calibrations between all 
of these samples do not match exactly, but they suffice to give the 
reader an idea of how our results compare with past work. In par- 
ticular, in comparison with previous work we note that while the 
2SLAQ data show an excess at 20 < g < 20.6, it generally shows 
a deficit for g > 20.6. The one exception is the faintest z < 2.2 
point from Koo & Kron (1988); however, that sample has a lower 
redshift limit of z ~ 0.9, whereas our sample extends to lower 
redshift. Overall, to the limit of our bright sample (g < 21.15), our 
agreement with previous work is well within the errors. Fainter than 
g — 21.15, if anything the 2SLAQ counts are deficient, but are still 
consistent considering the large coverage and spectroscopic com- 
pleteness corrections at these limits. Figure fTTI shows the cumula- 
tive 2SLAQ and 2QZ/6QZ quasar number counts. At the limit of 
the 2SLAQ survey, the cumulative counts compare well with the 
J = 22 cumulative counts (86.3 ± 16.5) from Zitelli et al. (1992). 
The slope of the cumulative counts are given as 3-bin averages by 
the dashed lines and the numbers at the bottom of the plot. The 
brightest 2SLAQ points are unreliable as 2SLAQ does not include 
quasars brighter than g = 18. The cumulative 2QZ/6QZ number 
counts gives a better idea of the slope at the bright end. Table [2] 
shows a comparison of the cumulative number counts predicted by 
the Boyle et al. (2000), Croom et al. (2004) and 2SLAQ best fit 
maximum likelihood parameterizations (assuming a double power- 
law and luminosity evolution characterized by a 2nd order polyno- 
mial) for g > 16.0 and 0.3 < z < 2.2. 



5.3 Luminosity function 

Figure fT2l shows two determinations of luminosity function derived 
from our sample. We first use the Page & Carrera (2000) implemen- 
tation of the 1/V method (Schmidt 1968; Avni & Bahcall 1980), 
which is shown by the points with error bars. This implementa- 
tion corrects for incompleteness at both the bright and faint limits 
of the survey. These incomplete bins (those not filled in Figure |8j 
are shown as open rather than closed points to indicate that they 
have been corrected for partial coverage of the bin. However, we 
note that the Page & Carrera (2000) correction for incomplete bins 
is not fully accurate since the (relatively large) z — M g bins are 
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Figure 11. Cumulative 2SLAQ (black circles) and 2QZ (red triangles) as a 
function of g. The numbers at the bottom indicate the slope of the 3-bin least 
square fits shown by the series of dashed lines. Note that the bright limit of 
the 2SLAQ data is g = 18 which causes a deficiency in the cumulative 
number counts at the bright end; at the faint end this lack of bright quasars 
makes little difference. While there is no strong characteristic break, the 
number counts clearly flatten with fainter magnitude. The blue pentagon 
shows the cumulative J-band number counts from Zitelli et al. (2002). 



not uniformly sampled, see Figure|8| The redshifts are the same as 
those in Figure 20 of Croom et al. (2004) for ease of comparison. 
The size of the redshifts bins is Az = 0.283 and the z = 1.39 data 
are repeated as grey lines in each panel. 

We next give the luminosity function as derived from a max- 
imum likelihood analysis; these are plotted as dashed/solid lines, 
the dashed part indicating extrapolation beyond the data used for 
the fit. The cyan lines show the best fit double power-law model 
(see below) from row 1 in Table 6 of Croom et al. (2004), which 
provide a poor fit at the faint end. The yellow lines show a similar 
model from row 1 in Table 3 of Boyle et al. (2000) (corrected to 
our cosmology), which has a steeper faint-end slope. Our own fit is 
shown in red and was derived as described below. 

We have assumed a luminosity function in the standard form 
of a double-power law (Peterson 1997; Croom et al. 2004) 4 



*( L A ■ 



or 



®(M g ,z) = —pjT——; 



10 0.4(a+X)(Ar s -AfJ) + 10 0.4(/3+l)(M 3 -A/ 9 *) ' 



(5) 



(6) 



We assume that the evolution with redshift is characterized by pure 
luminosity evolution (individual quasars getting fainter from z = 

4 We remind the reader of the well-known sign error in Boyle et al. (2000) 
whereby (in the convention used herein) the first equation in Section 3.2.2 
of Boyle et al. (2000) should have negative signs on a and /3 and the entries 
for a and fi in Tables 2 and 3 should be multiplied by —1. In addition, equa- 
tion 10 in Croom et al. (2004) and the equivalent equation in Section 3.2.2 
of Boyle et al. (2000) are missing a 1/L* factor in the numerator. 
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Figure 12. Faint quasar luminosity function for 2SLAQ quasars. Redshift bins are the same as for Croom et al. (2004). Open points are incomplete bins (see 
Fig.|8j. Corrections for photometric, coverage and spectroscopic incompleteness have been applied. Only the quasars within the 41 new plate centers with 
0.4 < z < 2.1 are included. The z = 1.39 data are repeated as grey curves in each panel for comparison. Best fit models (luminosity evolution characterized 
by a 2nd order polynomial in redshift) are shown from Boyle et al. (2000) [yellow], Croom et al. 2004 [cyan] and this work (red) as solid and dashed lines. The 
dashed parts of the best fit lines indicate where the fits have been extrapolated beyond the data. The faint quasar data suggest a steeper slope than the Croom 
et al. (2004) models. 



2 to today), with the dependence of the characteristic luminosity 
described by a 2nd-order polynomial in redshift as in Croom et al. 
(2004) where 

M* g (z) = M;(0) - 2.5(ftiz + k 2 z 2 ). (7) 

Note that this form assumes symmetric redshift evolution about the 
peak. This assumption is appropriate for UVX samples such as this 
one, but will break down for samples that extend to higher redshifts 
(e.g. Richards et al. 2005). 

We compute the maximum likelihood solution via Powell's 
method (Press et al. 1992) using the form given by Fan et al. 
(2001b). We first attempt to determine the best fit parameters by 
allowing all of the parameters to vary. The resulting parameters are 
given in the last row of Table 2 and the fit is given by the red line 
in Figure fT2l Due to the large incompleteness in our last magni- 
tude bin, we have performed these fits to a limiting magnitude of 
g < 21.65 rather than g < 21.85. The errors on the parameters are 
a a = 0.2, a = 0.03, a M * = 0.09, a kl = 0.02, a k2 = 0.01. 

Since there are relatively few bright quasars in our sample to 
tie down the bright end slope, we have also attempted to fix all of 
the parameters except for the faint end slope (/3) and the normaliza- 
tion to those found by Croom et al. (2004), specifically a — —3.31, 
M* = -21.61 + 51og/i 7u , fcl = 1-39, and k 2 = -0.29. 



The resulting faint end slope is then f3 = —1.45 ± 0.03 (with 
$* = 1.83 x 10~ 6 /i? Mpc~ 3 mag- 1 ). For both of these fits, a 
X 2 comparison of this model to the 1/V a data is formally rejected; 
see Table 2. We also note that there is apparently substantial co- 
variance between the parameters. For example, there is a signifi- 
cant difference in the faint end slopes of the Boyle et al. (2000) 
and Croom et al. (2004) analysis (as shown by the cyan and yel- 
low lines in Fie. ll2i . yet there is only a 1 per cent difference in the 
total expected counts to the limiting magnitude of the 2QZ survey 
(bj = 20.85). To the fainter limit of our survey, we find that the 
final 2QZ parameterization (Croom et al. 2004) significantly under- 
predicts (by 32 per cent) the total number of quasars to g < 21.65, 
while the Boyle et al. (2000) parameters yield a much better fit to 
the 2SLAQ data (see Fig. ll2l and Table[2j. The deviation from the 
best fit 2QZ model can be seen better in the left-hand panel of Fig- 
ure !13l where we have normalized our derived values by the best fit 
polynomial evolution model from Croom et al. (2004). The right 
hand panel is similar except that the data have all been normalized 
to our z = 1 .39 model in order to better show the redshift evolution 
of the quasars. All of the above suggests that the adopted parame- 
terization is not the optimal one; however, it still has considerable 
utility in terms of predicting counts of faint quasars and as an input 
for theoretical models. 



© 2005 RAS, MNRAS 000,[T]-?? 



12 



G. T. Richards et al. 



2.5 



"oj 2 



a 




T) 









2 




XI 




o 




rv 


1 


0) 




w 




^3 




CD 












BO 








>& 




u 




> 




— 1 




jO 




"o 





K 




-22 



-24 



-26 




Figure 13. Ratio of luminosity function from Fig. ll2l to the polynomial evolution models from Croom et al. (2004) . Colours and points are as in Figure 20 
of Croom et al. (2004) , specifically blue circles, green triangles, yellow squares, red circles, blue triangles, and green squares, corresponding to redshifts 0.54, 
0.825, 1.11, 1.39, 1.67, and 1.955, respectively. Note the deviation from unity at the faint end in each redshift bin. Right: The ratio of the luminosity function 
to our z = 1.39 maximum likelihood model. 



Table 3. Summary of maximum likelihood fits for the parameterization adopted for the 2QZ analysis (double power-law with luminosity evolution param- 
eterized as a 2nd order polynomial in redshift) and our adopted cosmology (§ 1). The redshift limits is 0.4 < z < 2.1, and objects must be brighter then 
M < —22.5. Nq indicates the number of quasars per square degree expected for 18.0 < g < 21.65. 



Sample 


a 


& 


M* 


fel 


k 2 


$* 


Nq 


x 2 


V 




Boyle et al. (2000) 


-3.41 


-1.58 


-21.92 


1.36 


-0.27 


9.88e-7 


66.8 








Croom et al. (2004) 


-3.31 


-1.09 


-21.61 


1.39 


-0.29 


1.67e-6 


54.4 








2SLAQ + Croom et al. (2004) 


-3.31 


-1.45 


-21.61 


1.39 


-0.29 


1.83e-6 


83.8 


161.5 


55 


2.1e-12 


2SLAQ only 


-3.28 


-1.78 


-22.68 


1.37 


-0.32 


5.96e-7 


79.8 


149.0 


51 


1.5e-ll 



We have also attempted to use the parameterizations of the lu- 
minosity function that were used by Wolf et al. (2003) since our 
data, like that of COMBO- 17, appears to show less of a break in 
the luminosity function than previous work. The best fit forms and 
parameters from Wolf et al. (2003) match the 2SLAQ data over 
a limited range in redshift and absolute magnitude, but these fits 
do not agree with the 2SLAQ data at the bright end and for lower 
redshifts. We were also unable to derive better fits to the 2SLAQ 
data using such parameterizations, likely because of the lack of dy- 
namic range at the bright end of the distribution. However, it is 
clear that other parameterizations, like those adopted by COMBO- 
17, are worth pursuing. 



5.4 X-ray comparison 

We can test the robustness of the faint end of the 2SLAQ luminosity 
function by comparing to quasar luminosity functions derived from 
X-ray selected samples which are thought to suffer less incomplete- 
ness as a result of the dust-penetrating nature of X-ray photons. In 
Figure [141 we compare the z — 0.825 and z — 1.67 redshift bins 
from Figure fl2l to Croom et al. (2004) and two quasar luminos- 



ity functions derived from hard X-ray surveys (Ueda et al. 2003; 
Barger et al. 2005). 

In these comparisons, we have converted between M g and 
log(L[2 — 10keV]) as follows. First we take our K-corrected M g 
and convert it to (rest-frame) log (L g [erg s _1 cm -2 Hz -1 ]) as pre- 
scribed by the definition of an AB (Oke & Gunn 1983) magnitude 
and an absolute magnitude (with an assumed distance of 10 pc) 

log(L g ) = -0.4(M 9 -5+48.6)+log(47r)+21og(3.086xl0 18 ).(8) 

Next we assume a power-law spectral index of ct v — —0.5 to con- 
vert from rest-frame g (4669 A) to rest-frame 2500A according to 

log(L 25 oo) = log(L 9 ) - 0.51og(4669/2500). (9) 

We then extrapolate to \og(L,2kev) assuming a luminosity depen- 
dent 2500A to 2keV slope, a ox , (Vignali, Brandt & Schneider 
2003): 

a ox = -0.11 log(i 2 50o) + 1.85 (10) 
and 

log(L 2 feev) = log(Z/2500) + Ctox log ( _H*SK ) . (H) 

V ^2500 / 
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Figure 14. Comparison of the 2SLAQ optical quasar luminosity function 
with X-ray quasar luminosity functions from the literature. Shown are the 
z = 0.825 (triangles) and z = 1.67 (squares) QLF from 2SLAQ and the 
best fit models for those redshifts from three other papers. The models are 
from Croom et al. (2004) (blue solid line), Ueda et al. (2003) (dashed black 
line) and Barger et al. (2005) (red long dashed line). 

(Using the revised a ox — L2500 relationship from Strateva et al. 
2005 yields somewhat fainter X-ray luminosities [~ 0.15 dex] at 
the bright end of our sample and would slightly flatten the X-ray 
QLFs in Fig. 1141 ) Finally, we compute a 2-10 keV luminosity by 
integrating over the 2-10 keV range assuming a photon index of 
T = 1.9 (Qx = -0.9). For comparison with Ueda et al. (2003) 
we further correct for the fraction of X-ray type II AGN accord- 
ing to their Figure 8 and an optical type II fraction of 0.5, which 
is roughly consistent with their Figure 9. In our comparison with 
the broad-line AGN luminosity function of Barger et al. (2005), we 
have treated their parameterization as if it were for a 2-10 keV lu- 
minosity rather than a 2-8 keV luminosity (since we are primarily 
concerned with the comparisons of the QLF slope), and we have 
applied a correction factor of 0.5 in the overall normalization. Fur- 
thermore, our comparison with Barger et al. (2005) differs some- 
what from their comparison with Croom et al. (2004) in that Barger 
et al. (2005) converted the optical and X-ray luminosities to bolo- 
metric luminosities in a manner which assumes a constant a ox , 
whereas we assumed the luminosity-dependent a ox given above. 
For both comparisons with X-ray QLFs we have converted the pa- 
rameterizations to the cosmology adopted herein. 

For the sake of facilitating the comparison of optical QLFs 
to X-ray QLFs, we note that, in the syntax used by Ueda et al. 
(2003) in their Equation 6 (and similar notation used by Barger et 
al. 2005, Equation 1) <f> M = A/2.5, a = -(71 + 1), and f3 = 
— (72 + 1), where <f>M, a, and (3 are defined as in Peterson (1997), 
Equation 1 1.33 (and similarly by Croom et al. 2004, Equation 1 1). 

In each case, the X-ray luminosity functions show less curva- 
ture in the faintest 2SLAQ bins than does the best fit model from 
2QZ. This comparison is not meant to be strictly quantitative since 
X-ray selected samples are more sensitive to obscured quasars and 
the conversion between Mb and Lx involves a number of tenuous 



assumptions. However, these comparisons confirm that the steeper 
2SLAQ faint end slope, while based on large correction factors, is 
quite reasonable. In particular, the agreement with the results of 
Barger et al. (2005) is remarkable. 



6 DISCUSSION AND CONCLUSIONS 

We have compiled a sample of 5645 quasars with 18.0 < g < 
21.85 and z < 3 using imaging data from the SDSS and the spec- 
tra from the 2dF facility at the AAT. We find a clear turnover in the 
optical number counts; a single power-law is not a good fit over the 
magnitude range sampled. For 20 < g < 20.6, the 2SLAQ num- 
ber counts show a slight (but statistically insignificant) excess over 
previous surveys, but the cumulative number counts are roughly 
consistent with the faintest surveys to 22nd magnitude. 

In terms of the luminosity function, we find good agreement 
with the 2QZ results from Croom et al. (2004) at the bright end, but 
the faint end 2SLAQ data require a steeper slope (higher density of 
quasars) than the 2QZ results from Croom et al. (2004). The previ- 
ous 2QZ results from Boyle et al. (2000) agree significantly better 
with 2SLAQ at the faint end. The lack of a well defined character- 
istic luminosity and covariance between the maximum likelihood 
parameters can explain the good bright-end agreement between the 
parameterizations studied and the faint-end disagreement between 
2SLAQ and the final 2QZ results Croom et al. (2004). Compar- 
ing to type I quasar luminosity functions derived from X-ray sam- 
ples suggests that the slope of the faint end of the 2SLAQ QLF is 
more accurate than the extrapolated faint end slope of Croom et al. 
(2004). 

An understanding of the quasar luminosity function is an im- 
portant ingredient for many different types of extragalactic inves- 
tigations. In particular, as has been stressed by those working with 
X-ray selected samples, investigations that depend on the optical 
QLF explicitly may need to be reconsidered as a result of recent 
revisions in the luminosity function of unobscured AGNs (not to 
mention obscured AGNs). Many investigations have an explicit de- 
pendence on the optical QLF, for example Bianchi, Cristiani & Kim 
(2001) in their analysis of the UV background; Hamilton, Caser- 
tano & Turnshek (2002) in their estimate of the quasar host galaxy 
luminosity function; Yu & Tremaine (2002) in their investigation 
of the growth of black holes; Croom et al. (2002) and Wake et al. 
(2004) regarding the clustering of AGN; Oguri (2003) in his deter- 
mination of the expected number of lensed quasars; Richards et al. 
(2004b) in their assessment of the lensing probability of the most 
luminous high-redshift quasars; and Fan et al. (2001a) in terms of 
the evolution of quasars from z = to 2 = 6. The QLF has taken 
on even greater importance in recent years with the realization that 
most massive galaxies host supermassive black holes, the correla- 
tion between black hole mass and stellar velocity dispersion (e.g. 
Magorrian et al. 1998; Ferrarese & Merritt 2000; Gebhardt et al. 
2000), and the possibility that feedback from quasars may play a 
role in the evolution of galaxies in general (e.g. Begelman 2004). 
In particular, models like those of Kauffmann & Haehnelt (2000), 
Wyithe & Loeb (2002), and others that attempt to explain the evo- 
lution of galaxies and quasars, rely on comparison with accurate 
observational determinations of the QLF. 

In fact, although the optical QLF presented herein is arguably 
the most robust determination to date for a large optically selected 
sample, for many applications an X-ray or far-IR QLF is more ap- 
propriate. That said, the luminosity function of optically selected 
quasars will remain an important tool for extragalactic astronomy. 
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The primary reason for this is the sheer size of the optical quasar 
sample (likely over 300,000 in the current SDSS imaging data 
alone; Richards et al. 2004a). While the deepest X-ray surveys may 
uncover thousands of AGNs per square degree, they do so over 
only a fraction of a square degree and the sum total area of the 
sky covered by both Chandra and XMM-Newton is unlikely to ever 
exceed even 1 per cent. IR surveys with Spitzer will cover a some- 
what larger area than X-ray surveys, but not at nearly the same 
space density as in the X-ray or with nearly the same total num- 
ber as in the shallower, but much wider optical surveys. Thus, this 
sample of faint quasars and the luminosity function derived from it 
will continue to provide important inputs to future extragalactic in- 
vestigations such as the Dark Energy Survey (DES) and the Large 
Synoptic Survey Telescope (LSST). 
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